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I. INTRODUCTION 



Spectroscopy of positronium (Ps) provides a sensitive test of bound state theory based on 
the Quantum Electrodynamics (QED). Because of the small mass of electron and positron, 
the effects of strong and weak interactions are negligible compared with the accuracy of 
present experiments. For this reason positronium represents a unique system which can, in 
principle, be described with very high precision by means of the QED only. Tests of the 
QED predictions are made possible by the very high experimental accuracy of positronium 
spectroscopy 0. 

The gross spectrum of positronium is well described by the Schrodinger equation with 
the Coulomb potential. Energy levels are 



where n is the principal quantum number. For the purpose of interpreting modern exper- 
iments the precision of Eq. ([[]) is insufficient. Corrections to the energy levels can in part 
be described by the Quantum Mechanics; however, for a complete description one has to 
resort to the Quantum Field Theory (QFT). Unfortunately, an application of the QFT to 
the bound states is difficult and special methods have to be devised 0-f|. 

Various approaches to bound state calculations have been reviewed e.g. in 0. Here 
we focus on a method close to the so-called Non-Relativistic Quantum Electrodynamics 
(NRQED) [J5J, which is an effective field theory based on the QED, for small energies and 
momenta. Eq. (|l|) implies that the characteristic velocity of the electron and positron in 
positronium is of the order of the fine structure constant a C 1. It is appropriate to apply 
a non-relativistic approximation to this system. 

Recently much progress has been achieved in the framework of non-relativistic effec- 
tive theories, mainly by employing dimensional regularization. It has been shown |7j that 
this regularization procedure permits an exact separation of effects arising at various char- 
acteristic energy scales. Using that method, which we will call dimensionally regularized 
NRQED (NRQED e ) the complete energy spectrum of Ps has been reproduced to order ma 5 
||. More recently, we have computed ma 6 corrections to the hyperfine splitting (HFS) of 
the Ps ground state ||, confirming one of previously obtained numerical results ||10|| . In the 
present paper we generalize that result to all 5* states, confirming [fll[], and compute also 
their spin independent shift at 0(ma 6 ) (obtained numerically in p"Tf). 

It is convenient to describe the energy of an nS state of Ps by dividing it up into the 
spin-averaged part and a part dependent on the total Ps spin (hyperfine splitting): 



where J is the total spin value of the Ps and s± are the spins of the electron and positron, 
respectively. One finds: 



E(J,n) = E mer (n) + s + s_E Ms (n), 



(2) 



1 (triplet state) : 
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J 



s + s 



+ 



J 



(singlet state) : 
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(3) 



s + s 
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Both the spin-averaged energy and the hyperfine splitting can be represented by series in 
powers and logarithms of the fine structure constant. In the lowest order E aver (n) = E{n) 
is given by Eq. ([1]), and E h{s (n) = Oijna 4 ). 

To order ma 5 the results for E avpr and Ey,h were found in |IT2Hlll. Those corrections 



have several sources: electron and positron charge radii and anomalous magnetic moments, 
vacuum polarization, two-photon exchange, two-photon annihilation and one-loop correction 
to the single-photon annihilation. 

Current accuracy of high precision experiments requires a complete calculation of the 
0{ma 6 ) corrections Ai? aver and AE h{s . 

The most precisely measured property of positronium is the ground state HFS, i.e. the 
energy difference between the two lowest states with total spin 1 and 0. Two best experi- 
mental values are 



Av = E(l 3 Si) - ^(l 1 ^) = 203 387.5(1.6) MHz, 



found in [|T^,|T§] and 



Av = 203 389.10(0.74) MHz, 



(4) 



(5) 



obtained in p7| . Another quantity of the experimental interest is the energy difference of 
2 3 Si and l 3 Si states P|: 



E(2 3 S 1 ) - E^Si) = 1 233 607 216.4(3.2) MHz. 



(6) 



The absolute accuracy of this measurement is clearly less impressive than that of the hy- 
perfine splitting. However, since met 6 = 18.658 MHz, a complete calculation of the energy 
levels at this order is warranted. 

At order ma 6 both A.E aver and AE^g can be written as 



AE — AE'rad + AE'annih + Ai? ra( j rec + AE rcc . 



(?) 



The logarithmic contributions at this order, 0(ma 6 In a), present in the annihilation Ai? ann ih 
and recoil AE rec corrections, were found first [|T9| , |20l . AE Ta a arises from the radiative correc- 
tions to the Breit potential at 0(a, a 2 ) pl|j22| . The three, two, and one-photon annihilation 

mi 



and 



contributions giving AE annih were found in |p3| , fl24 
annihilation radiative recoil contributions Ai? ra( j 
recoil corrections AE rcc were obtained in |T0|,|5|,|29[| for the HFS and in [1 I ] for E. 



respectively. The non- 
were calculated in |2~7LE3, while pure 



In this paper we present an analytic calculation of the recoil and radiative recoil cor- 
rections, AE'rad and AE'radrec, to energy levels of arbitrary nS positronium states. The 
rest of this paper is organized as follows: in Section |lj we discuss our method in general 
terms. Section [TTT] is devoted to the calculation of the HFS. Many technical details of this 



calculation are discussed there. In Section [TV] we present a calculation of the average energy 
E'aver. It is very similar to HFS, except that some additional operators contribute. Also the 
0(ma 6 ) radiative recoil corrections are discussed. Our results are summarized in Section 0, 
where also an overview of the theoretical and experimental situation is given and a complete 
analytic formula for the nS energy levels to order ma e is presented. 
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II. FRAMEWORK OF THE CALCULATION 



Before getting into details, let us describe the general framework of our calculation of 
the 0(ma 6 ) corrections to energy levels. 

First, we calculate an on-shell scattering amplitude for non-relativistic (v <C 1) particles 
to the needed order (the fact that v ~ a in Ps serves as a counting rule for contributions of 
various operators). In addition to the leading, single Coulomb exchange, this includes the 
relative 0(v 2 ) Breit corrections and also higher order 0{y 4 , av 3 ) terms. This non-relativistic 
amplitude is gauge invariant, and taken with a minus sign provides the potential for non- 
relativistic particles. 

Next, we use the ordinary quantum mechanical perturbation theory to find the correc- 
tions due to that potential; as unperturbed states we use the solutions of the Schrodinger 
equation with the Coulomb potential. We get the 0(ma 6 ) correction to energy levels as the 
sum of the first order correction due to 0(v 4 , av 3 ) perturbation and of the second order cor- 
rection due to the Breit Hamiltonian. Previously, this scheme was used for the calculation 
of the 0(ma 6 In a) corrections to the levels of S-states |30| and of the 0(ma 6 ) corrections 
to the levels of P-states [|31| ,|3^1. 

In the present calculation the result of the non-relativistic calculation is divergent. This is 
because also the short-distance ( "hard" ) corrections contribute. They arise from virtual mo- 
menta regions of the order of electron mass and cannot be obtained from the non-relativistic 
expansion. 

Our calculation is performed in the spirit of NRQED. We apply dimensional regulariza- 
tion, which offers technical advantages over more common techniques, based on the introduc- 
tion of an intermediate cut-off to separate the relativistic and non-relativistic momentum 
regions. Dimensional regularization makes the matching of the low-scale effective theory 
and the complete QED extremely simple. We find that in the sum of the short and long- 
distance contributions the singularities in the parameter]] e disappear and one arrives at a 
finite result. 

The spinor algebra in dimensional regularization requires some comments. In order to 
obtain the energy shift due to an operator Oi one has to calculate the trace of the form 
Tr ^^Oi^ , where \1/ is an appropriate wave function. The spinor parts of the relevant wave 
functions are 

* p = W 75 ' *° = W 7t 

for para and orthopositronium states, respectively. In the latter case, £ is the polarization 
vector (we average over its directions). The traces are calculated in a standard way in the 
/^-dimensional space. One encounters only even numbers of 75 matrices, and we treat them 
as anticommuting. 

Since the matrix elements involve the positronium wave function, it is easiest to calculate 
for the ground state (n = 1). However, once the corrections to the ground state have been 
found, there is a convenient way of finding them for excited states, with an arbitrary value 



Throughout the paper, we use the following notations: D = 4 — 2e and d = 3 — 2e. 
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of the principal quantum number n. Only the non-relativistic contributions have a non- 
trivial dependence on n. Their computation in dimensional regularization would be difficult. 
However, this task is simplified using other regularizations. Finally, we eliminate the cut-off 
dependence by requiring that for n = 1 the result matches the formula we found for n = 1. 
The freedom of choosing the regularization scheme simplifies considerably this part of the 
calculation. 



III. HFS OF THE POSITRONIUM GROUND STATE 



In this Section we present a calculation of the recoil corrections to the Ps ground state, 
A rec i? h f s . It is given as a sum of soft (non-relativistic) Eq. flllD and hard Eq. (ffOf) scales: 



331 In 2 17C(3) 



ArecEhfs = A nonre ii?hfs + Ahard-^hfs = rna — — In a + — — — 



6 



432 



4 



127T 2 



• (8) 



Those two groups of contributions are computed, respectively, in Sections |III A| and [Ti [ B 



Further, in Section 1H C| , we find a generalization of this result for radially excited states 
(arbitrary n): 



A rPC E h fJn) 



ma 



rv' 



-- (in- + ^(n) +je] + — ■ 
6 V n y ' ' J Yin In 2 432 



I 295 _ In 2 _ 17C(3) 5 



4 



127T 2 



(9) 



where $?(n) is the logarithmic derivative of the T-function and je — 0.577216 is the Euler 
constant. The n dependence of this result and its numerical value at n — 1 are in agreement 
with pi. 



A. Soft scale contributions 

We divide up the non-relativistic contributions to HFS into 6 parts: tree level Coulomb 
and magnetic photon exchanges, retardation, one-loop operators, and second iteration of 
Breit Hamiltonian which includes intermediate S and D wave states: 

A nonr el-Ehfs = Ac-^hfs + Am-Ems + A ret -Ehfs + Ai_i oop _Ehfs + As-Ehfs + AD-^hfs- (10) 

These partial results, given in Eqs. (^Jj, [!3|, |35|, |51], |B7D add up to 



A nonrel £ M s = |^ 2 (0) (~ - 41n(ma) + ^) . (11) 

In the remainder of this Section we discuss in detail how these contributions are calculated. 
According to standard procedure |3^] we identify the on-shell scattering amplitude, taken 



with the minus sign, with the matrix element of an interaction operator in the momentum 
representation. The soft scale contributions are calculated using the time-independent "old- 
fashioned" perturbation theory and the Coulomb gauge. Since this technique is not very 
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common, let us recall its basic ingredients. Exchange of a Coulomb or magnetic photon is 
described, respectively, by — ina/q 2 or — Anacti ® otj(o~ij — QiQj/q 2 )/^\q\- In the latter case, 
the denominator 2\q\ arises from the magnetic photon's phase space element. 

An intermediate state introduces the factor (E — E int + iO) -1 , where E int is the energy 
of the intermediate state and E is the total energy of the process. 

Dirac spinors are 

u(p) = J-^A + {p)w, (12) 
y u p + m 

where w denotes the four-spinor of a particle at rest; projectors on the positive and negative 
electron energy states are given by 

cxp + f3m \ r— 2 

1 ± , uj p = \Jp z + m z . 

u p J v 

In an expression for the potential the projector A_ contributes an additional minus sign. 

We begin with the contributions of the tree level effective operators, describing an ex- 
change of the Coulomb or magnetic quanta. The tree level operators, relevant for the 0(ma 6 ) 
calculation of the HFS, arise as 0(v 2 ) corrections to the Breit potential. 




1. Tree-level Coulomb photon exchange 



For the HFS we need the spin-dependent part of the 0(v 4 ) correction to the Coulomb 
exchange (see Eq. (|A2|) ) : 



V c (p',p) 



ixa [a-p, crp'] [cr'p, cr'p'\ 
16m 4 q 1 



(13) 



To calculate the spin part of the matrix element, we take the trace with <i-dimensional 
sigma-matrices and find (the factor 1/d in Eq. (|l4|) arises from the average over directions 
of the o-Ps polarization vector) 



1 d — 4 r 

— Tr (a i [*p,(ri/]a i [tr'pf,(T'p]) = 4—^- 



pV - (pp) 2 



-Tr{[(Tp,(Tp , \[a , p , ,(T'p\) =4 



P' 2 P 2 



{p'pf 



respectively for ortho and parapositronium. Using 

P' 2 p 2 -{p'pf = {pp) q 2 -{p'q){qp) 



(14) 
(15) 

(16) 



and noting that the average value of p'p in an S— state vanishes, we obtain the contribution 
of Vc(p,p') to the ground state HFS: 



A c E Ms = {V c (p',p)) 



5=1 = ^na_ I {p'q){qp) 
s=o dm 4 \ q 2 



(17) 



In Eq. fll7Q the matrix element is to be calculated over the ground state wave function in d 
dimensions: 

d d p d d p' 



(27r) d (2tt) £ 

Let us briefly explain how the integral in Eq. ([H]) is calculated. Although the integrand 
does not look complicated, the difficulty is that the exact form of the wave function ip(r) in 
d dimensions is not known. Fortunately, it turns out to be unnecessary. 

There are two alternative ways to calculate this integral. One is to transform it to the 
coordinate space. A divergence arises at r = and in the final result is proportional to the 
ci-dimensional ^(O); the remaining, finite part can be easily calculated in d = 3.0 

In the alternative approach we use the fact that the wave function in Eq. (|T7|) satisfies 
the (i-dimensional Schrodinger equation, which in the momentum space readsj] 

Anam r d d k 6(h) . 
^ P) = W^EJ (2n) d (p-kr (18) 

Using this equation we rewrite the integral in Eq. (O) as 

/ (P'q){qp) \ I (47f am) 2 (p'q)(qp) \ 

\ 9 2 /pjy \ (P 2 - mE)(p - k) 2 q 2 (p' 2 - mE)(p' - k') 2 / k k > 



(19) 



where the integration over p, p', as well as k, k', in the last expression is understood. The 
integral over p and p' receives a divergent contribution only from the region where p and 
p' simultaneously become infinite. Therefore, a single subtraction is sufficient to make this 
integral finite. It is convenient to subtract from (|T9|) the following expression: 

/ {Anam) 2 {p'q){qp) \ 



\ (p 2 — mE) 2 q 2 (p' 2 — mE) 2 / uu' 

After the subtraction is done, two nice features emerge. In Eq. (B^) the integration over 
k, k' factorizes and leads to ip 2 {0) times a two-loop integral, which can be easily calculated 
for arbitrary d. On the other hand, the difference between the last integral in Eq. (|19|) and 
the integral in Eq. ( p0| ) is finite and can be calculated for d = 3 using the explicit form of 
the wave function, 



I iram 2m 2 a 2 ma 



2 



We note that the counterterm (50) is constructed in such a way that the above mentioned 



difference vanishes for the ground state. This can be easily seen by integrating over k, k' in 



2 In general, also the derivative of the wave function at the origin, dip(r)/dr at r = 0, can appear 
in the divergent part of the integral. However, the Schrodinger equation relates it to i/>(0). 

3 This equation corresponds to a summation of an infinite number of ladder diagrams in the 
Coulomb gauge. For consistency it is essential to use here dimensional regularization in the same 
way as in the other loop integrations. 
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Eq. fl20D and using the fact that the p, p'-dependent terms in the denominator of Eq. (p0|) co- 
incide (up to a normalization factor) with the three-dimensional ground state wave functions 
in the momentum representation. 

Both methods described above lead to the same result. For d = 3 — 2e we obtain:^ 

AcEms = 2^ 2(o) 6 - 4 Hma) - 1) ' (2i) 

where ip(0) is the value of the d- dimensional ground state wave function at the origin. 



2. Tree-level exchange of a magnetic photon 



We now consider the correction caused by the tree level exchange of a magnetic photon, 
Fig. |l|(b). We neglect the energy dependence in the photon propagator; it will be restored 
in the following Section, where we discuss retardation effects. The relevant potential is 
obtained from Eq. 



V M (p',p) 



na [tr'q, a' 1 } 
16m 4 q 2 



P + [trq,(r\ [p +p 



J2 



Contribution of this interaction to the ground state HFS is 



A A f£hfs = (V M (p',p)) 



5=1 
5=0 



( d — 1 hoc 

d m 4 



12 



P +P + 



(P 2 -P' 2 f 



2q 2 



(22) 



(23) 



In d = 3 this matrix element is linearly divergent. To demonstrate how we treat linear 
divergences let us consider the p 2 term on the RHS of the above equation: 



P~ 



d d p 

w 



d p 2 (p(p) = mip(0) 



d d p 
(27p 



E(j){p) + 



d d k Ana 
{2n) d {p-k) 2 



<P(k) . (24) 



Shifting the integration variable p — > p + k we find that the p-integral in the last term is 
scale-less. In dimensional regularization such integrals vanish. The first term in Eq. fl24|) is 
finite in three dimensions. We obtain 



P~ 



m^ 2 (0). 



(25) 



Applying a similar procedure to the last term in Eq. ( |2"3"P we find the contribution of Vm(p', p) 
to the ground state HFS: 



M^Ms 



m 4 



mVf(0)-4 



d-1 l {p'q){qp) 
d \ a 2 



7ia 
3m 2 



^(0) - -41n(ma) + - . (26) 



4 We neglect factors T 2 (l + e) and (Att/j, 2 ) 26 which do not contribute to the final, finite result. 
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3. Retardation effects 



Let us now consider the retardation effects, which mean that the magnetic photon emitted 
by the electron propagates for a finite amount of time before being absorbed by the positron. 
During this time, the electron and positron can interact by several Coulomb exchanges 
(Fig. |l|(c,d,e)). To calculate the influence of these effects on the HFS, it is sufficient to take 
the spin-dependent parts of the current j(p',p) = u + (p')au(p) in the leading nonrelativistic 
approximation: 

The scattering operator describing the retardation effects is nonlocal both in space and time: 

k k 

a f d d k , Atr'k, cr-]47r <>y ~ TF [<rfc,cr,-] . , TT , . 

Here we assume that the magnetic photon with the momentum k is emitted by the electron 
at a point r e and absorbed by the positron at a point r p . Between those moments, the 
evolution of the system "positronium + photon" is governed by the propagator 4Tr(5ij — 
kikj/k 2 )/(2k)(k + H — E)~ x , H being the Hamiltonian of the nonrelativistic positronium 
slowly moving due to recoil. In the region of interest k ^> E and one can expand the 
amplitude ( p8|) over the powers of (H — E)/k ~ a. The zeroth term of this expansion is the 
spin-dependent part of the Breit potential, 

-AS(<,) = ~ l "' q ' a "^ q - a ' ] . (29) 

We need the second order term: 

f d d k 4:71 [er'fc, a h ] rTT , , NirTT ,., .Aak.a 1 ] , . 

V vet = aj ^^L_^ [F) e Xp (-^)][i?,ex P ^r e )]i 1 ^+H.c. (30) 

Only kinetic part of the Hamiltonian, 

H kin = + Ei. (31) 
km 2m 2m' y J 

has to be retained in the commutators. We find 

d d k 4vr [cr'k,a H ] n ^ , OJ _ , f .^ f _ _ , wjJ2 , nu _ s [trk,<r 

(27r) d 2fc 4 8 



Ket = -ay ^^4 L -^T J ( fe + 2fc P P ) ex P (»*(»"e - r p )) (fc 2 + 2fc Pe )^^ + H.c. 

(32) 



Transforming back to the relative coordinate r = r e — r p and the relative momentum 
p = p e = —p p , we get for the ground state HFS: 

^ net ( m 2 a 2 , 9 . , d—ll (p'q) (qp) \ \ rra 3 , 9 . . / 1 , . , 1 \ 
A re t-Ehfs - ^ *»(0) - 4— (^«) j = ^ 2 (0) (- - 4 ln(ma) - -) . 

(33) 
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4- One-loop operators 



Now we turn to the operators generated by one-loop diagrams. For the HFS the only 
contribution comes from the graph in Fig. |l|(f), which describes the mixed Coulomb- magnetic 
exchange with a transition of one of the particles to a negative energy state. In other words, 
this corresponds to a creation of an additional electron-positron pair by the electric or 
magnetic field of the electron or positron. 

Using Feynman rules for the time-independent perturbation theory, given at the begin- 
ning of this Section, we derive the corresponding potential: 

It induces the following correction to the ground state HFS (<i-dimensional integration 
over k is implicitly assumed below) 

J- la 2 / Antf - k) 4n(k - p) \ 
^i-ioop^hfs- d m3 \( p ,_ k) 2 {k-pY I 

- 1 a 2 / p' 2 +p 2 An 47r Anp' 1 Anp 1 \ 



d m 3 \ 2 (p> - k) 2 (k -p) 2 {p' - k) 2 (k - p) 2 / 
l7Ta ij 2 (0) (-- A\n(ma) - -) . (35) 



3m 2 V e 3 



5. Breit Hamiltonian 



To complete the calculation of the soft scale contributions to the HFS we have to consider 
the second iteration of the Breit Hamiltonian. It is obtained by including the effects of tree 
level Coulomb and magnetic photon exchang well correction to the kinetic energy. 
Using Eqs. and (H3) we find 



U(p',p) 



P 1 



Am 3 



(2n) d 5(p' -p) + — + 



ira Ana (p'q)(qp) — (p'p)q 2 net [trq, a 1 } [tr'q, a 



m 2 m? 



Am 2 



(36) 



In the position representation this Hamiltonian becomes 

d — 1 f p 2 „ , , 1 dira „ , , 1 



U(r,p) 
where 



P 



+ 



Am 3 Am m 



-iC(r)} + —5(r) 



16m 2 



C(r) = 



aT(d/2 - V 

n d/2-l r d-2 



[<TVXKV,a'*],C(r) , (37) 



(38) 



is the d-dimensional Coulomb potential. 
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6. Second iteration of the Breit Hamiltonian: S-wave 



We consider first the contribution of the intermediate S— states. The S"-wave part of the 
Breit Hamiltonian (571) reads 



U s (r,p) 



P 1 



+ 



4m 3 4m 



d — 1 { p 2 .1 duct . . net , 
, C (r) \ + - r Hr)-— | 



m 



a 4 ,a,]K,a']5(r). (39) 



It is convenient to divide up the calculation of the Us contribution to the HFS into two 
parts and consider the first and the last two terms in Eq. (p9|) separately. We begin with 
the latter, which we denote by A S iE h f s : 



Sl-C'hfs 



, {d - l)(3d - 2) ^na j ^ ^ ^ 



d 2 



, |m(r'))(m(r)| 
-E — £L 



5(r) 



8 



(d-l)(3d-2) /vra 



d 2 



, KQ)P 

E — E r , 



(40) 



In three dimensions the last sum is ill-defined due to ultraviolet divergences in the zeroth and 
first terms of its expansion in a. We denote these singular terms by G o (0, 0) and Gi(0,0), 
respectively, and obtain 



Aq-iE- 



Sl-^hfs 



,(d-l)(3d-2) ^\% 2(0) 



d? 



77r 



3m 2 a 



8tt 



+ G (0,0) + Gi(0,0) 



(41) 



Go(0,0) and Gi(0,0) are calculated in d dimensions, 
d d p 



G„(0,0) = - 
G x (0,0) = - 
and one finds 



m 



m 2 a 



{2n) d p 2 -mE 8tt ' 

d d p' r d d p m Ana m 



m 2 a f\ 



(2n) d J (2n) d p' 2 -mE q 2 p 2 - mE 16vr Ve 



AsiE h fs 



7na 3 (\ A1 , , 115 

- mi ' - 41n(ma) H 

,e 21 



9m 2 



•^ 2 (o) 



-41n(ma) + 2) , (42) 



(43) 



The contribution of the first two terms in Eq. (|39|) is calculated in the following way. We 
first write them as 



+ 



d-l \p 2 



Am 3 Am 



m 



C(r) 



1 

4m 



# 2 - d {H, C{r)} + (2d - l)C 2 {r) 



(44) 



where H = p 2 /m + C(r) is the leading order Hamiltonian. Correction to the HFS induced 
by Eq. (E3) reads 



d — 1 7ia 



d m 3 

We introduced here the reduced Green function 



d {H, C(r)} - (2d - l)C 2 (r) G(r, r')5(r') + H.c. 



(45) 
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G(ry) = E >|M|M (46 ) 

which satisfies the equation (H — E)G(r, r') = ip{r)ip{r') — 5{r — r'). Using obvious short- 
hand notations one can rewrite Eq. ([45]) as follows: 



A S 2^hfs = —- — n ( 6aE-G + 3-^(r)^(r') + 5— (G - G c 
Jm 1 \ \ r r r z 

3(d-l)(2d-l)^ 2/ ^ 



+ 2 " ^ -C 2 (r)G J *(0 + H.c. > , (47) 

We dropped massless tadpoles and separated the contribution of Go, which is the only one 
we have to calculate keeping d ^ 3. We find 

C 2 G 5(r') +H.c.) =— V 2 (0)Gi(0,0), 

-Gtf(r') + H.c.\ = -a<9 a (5(r)) = -3^ 2 (0), 



r 

o 
r 



-2£, 



- G )5(r') + H.c. ) = -4m^ 2 (0). (48) 



To obtain the last line we used the following equation: 

G{r, 0) - G (r, 0) = ~^^ r (hi(2 7 r) + lE - - + 7 r J , (49) 
where 7 = ma/2. ^From Eq. fl4"7| ) we now find 

A 52 £ hfs = f||V(0) (J - 4 ln(ma) + . (50) 

The sum of A^i^fs and As 2 -Ehfs gives the final result for the correction to the ground 
state HFS induced by the second iteration of the 5*-wave Breit Hamiltonian: 

A s ^ h f s = ^^ 2 (0) [- - 41n(ma) + —J . (51) 



7. Second iteration of the Breit Hamiltonian: D-wave 

Because of the last term in Eq. (|36D Breit Hamiltonian has non- vanishing matrix elements 
with |AL| = 2. In our case this causes virtual transitions from the triplet S'-state into in- 
states (transitions from the singlet state are forbidden by the total angular momentum 
conservation). Again, power counting shows that only the zeroth and the first order terms 
in the Green function expansion in a diverge in three dimensions. We first compute the 
remaining, higher order terms, which are finite for d = 3. 
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The sum of those higher order terms can be written as 

A^°£ hfs = (U D G CGCG U D ) , (52) 

where 

jj = a 3(crn)(cr'n) - era' 

Aw? r 3 

is the |AL| = 2 part of the Breit Hamiltonian in three dimensions, G and Go are defined in 
the previous Section, and C = —a/r is the Coulomb potential. 
The correction to the ground-state wave function, 

6 i/>(r)=G U D il>(r), (54) 

which appears in Eq. fl52]), satisfies an inhomogeneous Schrodinger equation: 

( E -^-)5o^(r) = U D ^(r). (55) 
y m J 

Solving this equation for 5oip(r) we obtain 

A h D -°E his = 85 S1 (£) (^G D (r, ri )^-) , (56) 
y 24/ \mr z mrf j 

where Go(r, r\) is the D-w&ve part of the Green function G, and the factor 85 si arises from 
Jcrcr' - 3(<xn)(<r'n)] 2 ) = (3 + 4a a' + (<t<t') 2 ) = 88 S1 . (57) 
To calculate the matrix element in Eq. fl5"6| ) we note that 

-L = hH D -H), (58) 
mr z 6 

where Hd(H) is the radial Hamiltonian for Z)(S")-states. Using equations of motion for both 
the Green function and the wave function in Eq. (|56|) one finds 



in agreement with [19 



To complete the calculation of the .D-wave contribution we have to consider the zeroth 
and first order terms in the a expansion of the Green function, 

A D0 £ h fs = (U d GqU d ) , (60) 
A D1 £ hfs = {U D G 1 U D ) . (61) 

The perturbation Uo{p',p) is extracted from Breit Hamiltonian, Eq. (1361), and reads 



13 



U»V,p) = {^Mlfl - MpL^] . (62) 
Am 1 \ a q J 

The average is taken over the (i-dimensional wave function. Calculating the trace using the 
triplet wave function we obtain 

// fogKofl [aq^\aW^ \ ( \a k MK^ jcrq,CT k ]lcT>q,a> k ] \\ 
\{ d q» ){ d <? 2 J I 

= {BijWB^q)) , B, tJ (q) ^{m-h\ (63) 



Therefore 



aid - 2) 2 

aid - 2) 2 I 4vr \ 

^DiEu* = ~ \ m J {B^p 1 - Ifigffljp-^gWBvik-p)} , (65) 



where 



g(k) = -^- 2 , (66) 

and (i-dimensional integrations over k in fl6"4|) and over fc, A;' in (|65l ) are understood. Some 
details of the integrations in Eqs. (64,65) are given in Appendix [B]. Adding the higher-order 
effects found in Eq. ([59]) we obtain the complete .D-wave contributions to HFS 

5vra 3 , 2 /l , x 19 



72,/ (0) r 4h(ma) "il- (67) 



B. Hard scale contribution 

Another contribution to the HFS arises from virtual momenta scales of the order of the 
electron mass. It can be calculated by considering the on-shell e + e~ scattering amplitude 
with an exchange of three photons in the t-channel (see Fig. 0) exactly at the threshold, 
i.e. for zero relative velocity of the incoming electron and positron, in dimensional regular- 
ization. The use of the dimensional regularization brings in essential simplifications, since 
almost any other regularization would bring in power-like divergences and hence require 
additional subtractions. This so-called hard scale contribution gives rise to four-fermion 
operators in the low-scale Lagrangian or, equivalently, to the S(r) terms in the effective 
quantum mechanical Hamiltonian. 

Technically, this calculation is similar to the derivation of the matching coefficient of 
the vector quark- ant iquark current in QCD and its NRQCD counterpart, described e.g. in 
f3~4] , |3"5]| . Here we outline the main steps of this calculation. 
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An arbitrary Feynman integral which contributes to the hard scale part of the calculation 
can be written as 



T( , fd%d\ 1 



where 



Si = k 1: S 2 = k 2 , S 3 = (ki - k 2 ) , S 4 = k 1 + 2pkx, 
S 5 = k\ + 2pk 2l S 6 = kf - 2pk Xl S 7 = k\ — 2pk 2 , 

S s = (h - k 2 f + 2p{h - k 2 ), S 9 = (h - k 2 f - 2p(h - k 2 ), (69) 

and ai, . . . , ag are integers. In practice we encounter diagrams with only at most 6 different 
propagators, so that at least 3 exponents are zero. Applying the integration by parts 
technique ]36| to an integral /({aj}), one obtains a set of relations among integrals with 



various values of indices {a^}. Using these relations one can express any /({a^}) in terms of 
a few master integrals. This is most easily done using symbolic manipulation programs. 
The result for the hard scale recoil corrections (Fig. to the HFS reads || 

A hard £ hfa = ^ 2 (0) [-\ + 41nm - + H_ 6 ln2 > ). (70) 



e 7r* 7r 



C. HFS for excited S states 

The result for the HFS of the ground state can be used to obtain the HFS for an arbitrary 
excited state. The non-trivial dependence on the principal quantum number n arises only 
from the soft scale contributions. Therefore, one has to repeat the quantum mechanical 
calculation of the non-relativistic part using any convenient regularization (we use a cut-off 
at 1/m <C ro <C 1/ma) and compare the result with the known formula for n = 1, Eq. ([|). 
One finds 



ma 6 |\ .. , 1 / a „ „ s \ 7 1 



n 3 



[div]--[\n- + V(n)+iE) + 



6 V n J \2n 2n 



,2 



(71) 



The quantity [div] in the above equation stands for the unknown and n-independent con- 
stant, easily determined by requiring that for n=l Eq. (|8|) is reproduced. We then obtain the 
final result for the recoil corrections to the HFS splitting for an arbitrary nS state, Eq. ([5|). 

IV. SPIN- AVERAGED ENERGY LEVELS 

To obtain 0(ma 6 ) corrections to the triplet and singlet energy levels separately, we have 
to calculate E wer {n) (cf. Eq. (||)). An appropriate formula for this calculation is 

p 3 -^triplet + -Esinglct d -Erf-plot + -^singlet 
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It is known [^,^,^] that the recoil corrections E aver do not contain In (a) at the order ma e . 
In dimensional regularization this means that the hard-scale and soft-scale contributions are 
separately finite. f\ 

Conceptually, determination of E aver is very similar to the calculation of the HFS dis- 
cussed above in detail. The only difference is that several new operators appear, which 
contribute to -E a ver but not to the HFS. 



A. The ground state average energy shift 

We begin with the correction to E aver induced by the relativistic corrections to the dis- 
persion law, Up = \Jp 2 + m 2 (Fig. 0(g)). Expanding u p in \p\/m, we obtain: 

u =m + ]—- it- + P 6 + (72) 
2m 8m 3 16m 5 

The last term induces a correction of the appropriate order: 

Adis P £ave r = (p 6 (2n) d 6^(p - p')) = -A^V(O). (73) 
8m° N 1 64 mr 

The 0(v 4 ) spin-independent part of the tree level Coulomb exchange amplitude 
(cf. Eq. (|A1[ ) and Fig. |I|(a); we neglect terms odd in p, whose average vanishes in an S- 
state), 

gives rise to the following correction: 

Wry 3 /I J\ 

A c £ av e r = ^^'(O) (- - 41n(ma) + -J . (75) 

Virtual transitions to negative energy states induced by the Coulomb exchanges, 
Fig. |l](h), generate an effective spin- independent operator 

V c 4r) = -^[p,C(r)f. (76) 

This operator describes the energy shift due to a creation of an additional e + e~ pair by the 
Coulomb field of either electron or positron. The resulting energy shift is 

A c _£ avcr = -^^ 2 (0) (i - 41n(ma)) . (77) 

The spin-independent part of the tree level magnetic exchange, Fig. |l](b), induces the 
following shift in the energy levels: 



Because of power-like singularities, this is not necessarily the case in other regularization schemes. 
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Am£ „ . - ( (p , + p2) _ ) . ^ (0) (i _ 4 Hma) + d . (78) 

To account for the retardation in the magnetic photon propagation, Fig. |I](c,d,e), we 
use the approach described in the HFS case. Our starting point is similar to Eq. (]28|), 
except that now the full expression for the currents must be used, rather than just their 
spin-dependent part. We obtain 

7ra; 3 /l \ 

A ret £ aver = ^^ 2 (0) - - 4 ln(ma) - 8 . (79) 

The next contribution comes from the exchange of two magnetic photons with creation 
of an additional e + e~ pair in the intermediate state, Fig. [l|(i). We find 

A M M-^aver = -^^(O) (- ~ 4 ln(ma) - 3 J . (80) 

We proceed further with the correction to E^cr, induced by the second iteration of the 
S'-wave Breit Hamiltonian. The calculation closely follows the HFS case. We arrive at the 
following result: 

ira 3 , , /M|2 /l 433 

The iteration of the D-wave part of the Breit Hamiltonian only influences the energy 
levels of the triplet state because of the total angular momentum (L + S) conservation. For 
this reason, to obtain the required correction to E aver it is sufficient to multiply Eq. (|67D by 
the factor d/(d+ 1). We find 



A^ aver = -^- 7 \m\ 2 ( - - 41n(ma) + ^ ) . (81) 



(2) 



57ra 3 , 9/ . (\ , , s 119\ . , 

^ 2 (0)(^--41n(ma)-— j. (82) 



aver ~ 96m 2 

It is easy to see that in 3 dimensions the spin-dependent operators do not contribute to 
E aver . However, since we work with divergent integrals and use dimensional regularization, 
this is no longer valid for d ^ 3. In this case an "anomalous" situation arises: spin- dependent 
operators provide contributions of the form (d — 3)/e to E aver , which are finite as e — > 0. 
Part of these contributions has already been accounted for in the corrections induced by the 
Breit Hamiltonian. The remaining contributions give 

15 HQ 3 

Aanom^aver = ■ (83) 

64 m l 

The hard scale contribution, Fig. 0, is calculated in the same way as for the HFS. One 
finds: 

Tra 3 ,2/ x A3 9C(3) 33 \ 

The sum of all contributions presented above provides the (D(ma 6 ) pure recoil correction 
to the ground state energy: 

ma 6 /901 11 3C(3)\ ma 6 . . . . 

A rec Saver = + ~^ + = 2.48688 . . . , 85 

8 \576 2vr 2 vr 2 J 8 V ; V ; 



in very good agreement with the numerical result of Eq. (20) in |Tl|], — in f-(2. 484(5)). 
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B. Energy levels for arbitrary n 



To generalize the result Eq. flS5|) for arbitrary n, we proceed according to the program 



outlined in Section [III Q We repeat the calculation of the soft-scale contributions to E aver 
for arbitrary n using a different regular ization scheme. Namely, we set d = 3 and cut off 
the divergent integrals over r from below at some r <C l/(ma). The transition to three 
dimensions simplifies the calculation. We find 

= ([*»] + _-_ + -). (86) 

The n-independent term in the above equation is regularization dependent. It cannot be 
determined by considering the soft-scale contributions alone. Nevertheless, by matching 
Eq. ( |S6"D to the shift for the ground state Eq. (|^), the "value" of the divergent constant 
[div] is completely determined. We obtain 

A „ , , ma 6 /83 11 3C(3) 69 8 2\ 

"-w ^-^(ra + ^ + ^ + e^-w + nj- (87) 

This is our main result for the recoil corrections to the energy levels of positronium. It 



agrees with the partially numerical result derived in 11 . 



V. RADIATIVE RECOIL CORRECTIONS 



So far in this paper we have been considering pure recoil effects. Another class of the 
0(ma 6 ) corrections to positronium energy levels and their HFS are the so-called radiative 
recoil corrections, where one of the three exchanged photons is created and absorbed by the 
same particle (see Fig. [3]). 

Our technique is very convenient for the calculation of these corrections. The key point 
is that at (D(ma 6 ) the radiative recoil corrections do not receive any contribution from 
the non-relativistic scales. Thus it is sufficient to calculate the diagrams shown in Fig. |3| 
(supplemented by the electric charge, electron wave function and mass renormalization) 
exactly at the threshold. For the same reason, the n-dependence of the radiative recoil 
corrections comes only from the 1/ra 3 behavior of the n^-wave function at the origin. Some 



details of this calculation are described in Section |IIIB| in the context of the HFS. We obtain 

>,6 / ■-> \ 70 i i 4 

+ - In 2 



A 



rad rec- C/ hfs 



ma 



A 



rad rcc -'-'aver 



rv 



ma 



rv' 



03) 
2tt 2 

9C(3) 
8tt 2 



79 41 
48 + 367 2 



+ 



97 



1025 



144 432tt 2 



(88) 
(89) 



respectively for corrections to the HFS and to the average energy, in full agreement with 
the analytic results of Ref. ||28|| . For completeness, we give here separately the contributions 
of electron vacuum polarization effects to radiative recoil corrections (they are 

included in (E 



A VaC pol rp 

^rad rec-^hfs 



ma 



n 3 9tt 



2 - 



A VaC pol 771 

^rad rec-^a 



ma 



Tl" 



1 

36 



27vr 2 



(90) 
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VI. SUMMARY AND CONCLUSIONS 



The main new results of the present paper are the analytic formulas (0) and ( jS7[) for the 
pure recoil 0(ma 6 ) corrections to the HFS and spin-averaged energy levels of positronium 
nS states. These recoil effects provide the last pieces needed to present complete analytical 
formulas for the total corrections to E aver and E^ B . We use the parameterization introduced 
in Eq. (|), 



E(J,n)=E a 



n 



- 5 jo E Ms (n) 



(91) 



and find 



E R 



n 



ma 
4n 2 
ma 5 



ma 
32n 3 



+ 



ma 
16n 3 



11 

An 



16 



14 



-6 In a - — In k (n, 0) + — In - + #(n) + j E 



a 



-In *(n) -7b + 

n 



2435 
" 432 



7 

n 



17 



3 

14103) 

4 7T 2 

69 



+ 



137 



(« 



68 

^2 



In 2 + 



37 , 7 
3 In 2 + — 

45 3n 

1421 



277T 2 



12n 2 16n 3 



(92) 



and 



-Ehfs(n) 



7 ma 4 

12 
4297 



+ 



mo 

7m 3 
/221 



ln2 + 



+ 



1 



3456 V144 2vr 2 



In 2 



n 3 
53 



5 

24 



a 



In - + m n ) +7s + 



?? 



1367 
648tt 2 



32tt : 



C(3) + f - - 



5?? 



85 
96^ 



(93) 



We have also recalculated the radiative recoil corrections, Eqs. (p8| , |89D , confirming recent 
result of Ref. p8 |. Let us make a technical remark. In dimensional regularization, used in 
this paper, the calculation of the radiative recoil corrections is particularly simple. Since 
there are no low-scale contributions to the radiative recoil corrections, it suffices to calculate 
corresponding Feynman graphs exactly at the threshold. No matching or subtractions are 
required. 

Formulas (^,^), together with P state energy levels given in Appendix |C], can be used 
to compute quantities which can be directly confronted with experimental data. We use the 
following values for the Rydberg |^T[ and fine structure |42j constants: 



Rr 



ma 



3 289 841960.394(27) MHz, 



a 



1/137.035 999 59(51). 



(94) 



In addition to the full corrections O (ma G ) we include the leading logarithmic terms 
0(ma 7 \n 2 a) found in ||3| for HFS, and in for the spin-averaged energy levels: 



A LL E( J, n) = - 



499 
To" 



+ 7(1- AS. 



JO 



ma 7 In a 
32vm 3 



'10- 



(95) 



For the most precisely measured quantity, the Ps ground state HFS, we find 
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Au = 203 392.01(46) MHz. 



(96) 



The 0(ma s ) recoil corrections to this observable have been subject of some debate. In the 
literature three different results have been reported |5| JlO|j29[] ,p| Our result for this correction, 



Eq. (§), evaluates numerically to ma e (— |lna + 0.37632^). This is in excellent agreement 
with Ref. [ lOfl , where for the non-logarithmic part of the correction a number 0.3767(17) was 



obtained. The framework of our calculation is similar to Ref. [|I(J . However, in that study a 
different regularization method was used. The agreement of the results gives us confidence 
in their correctness. 

Comparing Eq. (p6|) with the experimental results, Eqs. (|]|jD, we observe a significant 
deviation of the order of 3 — 4 experimental errors. It is not very likely that the uncalcu- 
lated higher order effects alone can account for this discrepancy. The size of the 0(ma 6 ) 
corrections gives no indication of bad behavior of the perturbative expansion. On the other 
hand, the leading logarithmic term O (ma 7 In 2 a) is sizable. A calculation of the subleading 
terms at this order remains an important theoretical challenge. 

For another experimentally interesting quantity, the energy interval of the IS" — 2S tran- 
sition, we get 

E(2 3 S 1 ) - £(l 3 Si) = 1233 607222.18(58) MHz. (97) 

in fair agreement with the experimental result, Eq. (P). 

Other quantities, for which high precision measurements have been made or are being 



planned, have recently been reviewed in ||28|| . In Table [I] we update the theory predictions 
for those observables. Our predictions are in good agreement with |]28|| . We have been able 
to decrease the error bars by including the analytical results (|9^,^K|) and the value of the 
leading quadratic logarithms fl9"o|). 

Finally we would like to comment on our error estimates. The errors due to uncertainties 
in the fine structure constant and the electron mass are well below 0.1 MHz level. The dom- 
inant theoretical error source is the uncalculated remainder of the perturbation expansion. 
Although formally ma 7 ~ 0.1 MHz, the leading 0(ma 7 \n 2 a) terms contribute —0.92 MHz 



to the HFS . It remains very important to calculate the remaining, non- leading terms in 
0(ma 7 ). For the present analysis we assume that the leading logs (D(ma 7 In 2 a) dominate 
the higher order contributions and take half their size as the theoretical error estimate. 

The spectrum of the nS and nP positronium energy levels is now known analytically, 
including effects 0(ma 6 ). Our calculation for the nS levels was made possible by new 
theoretical tools which have their roots in the recent perturbative calculations in high-energy 
physics. We hope that these methods will find further applications. 

The agreement between theoretical predictions and experimental results in Ps spec- 
troscopy is impressive with a few exceptions. One can only hope to find something new 
and unexpected by trying to put these exceptions in line with the overall picture. We look 



6 After our HFS calculation was completed, we were informed about an independent numerical 
calculation of the recoil corrections |45[| . Although that study is still in progress, its preliminary 
results seem to agree with Ref. [jlOj and the present paper. 
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forward to future improved measurements of positronium energy levels and their confronta- 
tion with QED. 
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APPENDIX A: TREE LEVEL ELECTRON POSITRON POTENTIAL 

We present here formulas for the potential arising from a single Coulomb or magnetic 
photon exchange between an electron and a positron, valid to 0(v 4 ). The virtual annihilation 
is not taken into account here. We also drop those terms which annihilate the 5-state wave 
function. These formulas, valid in <i-dimensions, are useful in the derivations of HFS and 
spin-averaged energy levels. 

For a single Coulomb exchange between two particles of opposite charges, Fig. |l|(a), the 
minus on-shell scattering amplitude is 

- A c (p',p) = -p(p',p)p(-p', -p), (Al) 

where p and p' are spatial momenta of the incoming and outgoing electron; q = p' — p; and 
the charge density is p(p',p) = u + (p')u(p). In momentum representation a single Coulomb 
photon exchange gives rise to the potential 

Anaf q 2 5 (p 2 - p' 2 ) 2 + 6q 2 (p 2 + p' 2 ) + q 4 + [ap\ ap] [a'p\ cr'p] \ 
Uc{p ' p) = ~^~[ 1 -4^ + oW /• 

(A2) 



In the leading nonrelativistic approximation Eq. (|A2|) gives the Coulomb potential. 

Next we consider a magnetic photon exchange, Fig. [TJ(b). We neglect retardation effects. 
The scattering amplitude is 

- A M (p>,p) = ^j^p'^^-p', -p) - ^\ , (A3) 

where j(p',p) = u + (p')cx.u(p) is the matrix element of the current. The resulting potential 
is 
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Um(p, p') 



ixa I 4 
m 2 q 2 \q 2 



na 



i l\2 2/2 

PP ) - P P 



2m 4 q 2 



(P 2 +P' 2 ) [{pp'f ~ PV 2 ] - i[<T<Z,^][<x"?X 
-^-=^ ([<rg, a 4 ] [<r'P, a'"] + [<tP, a 1 } [<r'q, o*\ ' 



(A4) 



These formulas are valid in the center of mass frame. We use p and p' to denote incoming 
and outgoing electron momenta, and q = p' — p, P = p' + p. The primed cr-matrices act 
on the positron spinor. 



APPENDIX B: USEFUL INTEGRALS 

In this Appendix we present various integrals which were useful in the calculations pre- 
sented in this paper. The following formulas have been used throughout the paper, especially 
for the tree level diagrams: 

2 2 

ma 



P 2 ) = -^^(o), 



4\ ™2^,2 



p \ ma ( 2 



\q 2 1 16 
' P 2 P >2 \ m 2 a 2 2 , (\ T 



\ q 2 I 4 
[qp){qp')\ m 2 a 2 2 /l 



^ 2 (0) f--41n(ma) + -) 



ip 2 (0) --41mW)-l . (Bl) 



\ <? 2 / 8 ' Ve 

In the remainder of this Appendix we describe some details of the .D-wave contribution 
to the second iteration of the Breit Hamiltonian. First, we rewrite Eq. (|64T ) in the following 

way: 

(Bqtf - kMVB^k -p)) = (b^p* - k)g(k)47r kikj ~ 2k ^+P^\ 



k - p)< 

B tJ (q)p lPj - 2B tJ (p> - k)9(k) ^^- Pj - B^p' - k)g(k)-^-^PiPj ) • (B2) 



Here and below we use the Schrodinger equation in the form 



Att 

<f>(p) = g(p) {p _ k)2 m- (B3) 



Similarly, by rearranging terms in (^) we get 

/ 47T 

iB^p' - k')g{k') (fc/ _ fc)2 g(fe)%(fc - p) 
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/ An An An 

= (BvV - k)9(k)j^—y 2 p iPj - 2 Bij ( P > - kWjj^^k^jj^Pj 

+Bij (p> - k')9{k') J ^ wg {k) J ^ wPlP }l , (B4) 
Using the symmetry with respect to p <-> p', we rewrite the first term in Eq. 



In the same way, the second term in Eq. ( |B2"1) is transformed to 

/8vr 4vr(fc-p)p 4tt / ,x An(ki - p { ) \ 

Considering the divergent part of this expression we find 

(t M^rtr) = 1 ir V(0) |G ° (0 - 0) + Gi(M)l • (B7) 

The sum of ([54]) and (|65[) reads 
A m E m + A mBhf , = » [Go( o,0) + G l( 0,0)] - + MlWw) 



4m 4 d \ d 3 q 2 

47r /n, A n^ik-Pi 



-Pj 



An An 
-2S y (p - fc )<?(*/) (fc /_ fc)2 ^( fc ) (fc _ p)2 ^ 

47T 47T \ 

- k')g(k')- (fl ^ g (k)j j ^ Pi p j \ . (B8) 



Only two terms here contain the logarithmic divergence: 
An 2 a\d-2) 2 na 3 



(Ci(0,0)) =^^ 2 (0)Q-41n(ma)-^, (B9) 



m 4 <i 2 



i^r \-^f^) - ss^o) U " 41n(rott) " tJ ■ < B10 » 

All other terms are finite and we compute them in three dimensions. Here we list some 
useful integrals (x = fc/7 and a(x) = arctan(x)): 

/" d 3 p pi 7^(0)A; 



(2vr) 3 (fc-p) 2 



23 



d 3 p PiPj 
(2Tr) 3 (k-p) 



;<t>(p) 



7^(0) 
k 



1 / kjkj 5{j 

2 [~w ~ y 



x 2 + 3 ./• 



2 -3 



x(x 2 + 1) 



-am + 



.1 



x 2 + 1 



d 3 p d 3 F g{k')B v {k' - p) _ 7^(0) ffc^- 4/ 



(2tt) 3 (2tt) 3 (k - k') 2 



k V fc k 



a(x) 



1 + - "5- 

/ AX 



(Bll) 



Using these above formulas in Eq. (|B8|) we find the final result for A^q-Ems + A^E^. 



APPENDIX C: P STATE ENERGY LEVELS 

In this Appendix we present formulas for the energy levels of P states, to order Oima 
Correcting some minor misprints in Ref. [31,32] one finds 

^, m ® 2 m ot A /13 11 \ ma 5 /4 16 , , . A 

ma 6 / 69 559 169 20677 3 , 9C(3) 13 



n 3 V 512n 3 4800n 2 4800n 432000 80 160vr 2 128tt 2 

r-./ in\ ma2 ma;4 /5 11 \ ma 5 /5 16 , . 

E(n Pl) = " 4rf ~ 4^ U " 16^ " W U + T ln fc ° (n ' 

ma 6 / 69 77 25 1 , n C(3) 179 493 

+^r~ -t-^t + tttt o - — — + — In 2 



n 3 V 512n 3 320n 2 192n 48 32tt 2 3456tt 2 17280, 

ma 2 ma 4 f A 11 \ ma 5 /25 16 



E{H Po) = " 4^ " 4^ U " IfcJ " W ll8 + T ln ^ 1 

ma 6 / 69 119 1 923 1 „ 3 >/oN 203 



77" 



3 V 512n 3 240n 2 3n 4320 8 16vr 2 ^ v ' 576?r ; 



_.. n . ma 2 ma 4 /2 11 \ ma 5 / 7 16 . , , 1 
E(n°Pi) = -— - — ^ - - — - - — ? — + — In k (n, 1 



4n 2 4n 3 V3 l&n) 87m 3 \lt 
ma 6 f 163 23 69 1 \ 

V4320 + 12077 2 ~~ 512n 3 ~ 12nJ ' ^ ^ 



71 3 



For the numerical evaluations we use the following values of Bethe logarithms ln [k (n, I) / R c 



ln[fc (l,0)/i2oo] = 2.984 128 555 765 498, 
ln [jfe (2, 0)/Roo\ = 2.811 769 893 120 563, 

ln [jfeo(2, l)/Roo] = -0.030 016 708 630 213. (C2) 
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TABLES 



TABLE I. Theoretical predictions for experimentally relevant positronium transitions. 



Transition Theory [MHz] 
2 3 Si - l 3 5i 1233 607 222.18(58) 

l 3 Si-lLS 203 392.01(46) 

2 3 5i - 2 3 P 18498.25(8) 

2 3 5i - 2 3 Pi 13012.41(8) 

2 3 S"i - 2 3 P 2 8626.71(8) 

2 3 ^ - 2 1 Pi 11185.37(8) 

2 3 ^ - 2 1 < S 25424.67(6) 
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FIGURES 
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FIG. 1. Non-relativistic corrections to HFS and spin-averaged energy levels: (a,b) Coulomb 
and magnetic photon exchange; (c,d,e) retardation effects; (f) mixed Coulomb-magnetic exchange; 
(g) relativistic correction to the dispersion law; (h,i) double Coulomb and magnetic exchange. 
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FIG. 2. Feynman diagrams representing pure recoil corrections to positronium HFS and 
spin-averaged energy levels. Wiggly lines denote photons in Feynman gauge. 
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FIG. 3. Examples of radiative recoil corrections to positronium HFS and spin-averaged energy 
levels. 
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